An Improved Harmonic Map Ansatz 
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The rational map ansatz of Houghton et al U| is generahsed by aUowing the profile function, 
usually a function of r, to depend also on z and z. It is shown that, within this ansatz, the 
energies of the lowest _B = 2, 3, 4 field configurations of the SU{2) Skyrme model are closer to the 
corresponding values of the true solutions of the model than those obtained within the original 
rational map ansatz. In particular, we present plots of the profile functions which do exhibit their 
dependence on z and z. 

The obvious generalisation of the ansatz to higher SU{N) models involving the introduction of 
more projectors is briefly mentioned. 
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Introduction 

A few years ago Houghton et al [j] presented an 
ansatz, the so-called rational map ansatz, in order to ap- 
proximate multiskyrmion solutions of the SU{2) Skyrme 
model by field configurations in which the angular de- 
pendence was determined by a rational map and the ra- 
dial dependence was determined by a numerical solution 
of a nonlinear ordinary differential equation. This last 
equation involved the so-called profile function (i.e. the 
generalisation of the one skyrmion profile function) and 
its shape had to be determined numerically. 

This ansatz has been a great step forward since it gave 
very good approximations to the solutions of the full 
equations which up to then could only be determined 
numerically (and these simulations involved hours, days 
or weeks of CPU time). In fact, these approximations 
were so good that the values of the energies were only, 
at most, a few % up on the true value and it was prac- 
tically impossible to distinguish the energy density plots 
obtained with the use of the ansatz from the exact ones 
obtained numerically. 

At the same time the ansatz clarified the situation 
for higher baryon number states; the energy density had 
maxima at several points lying on a shell whose radius 
grew with the baryon number. It has also lead to some 
generalisations - first, the generalisation to the SU{N) 
Skyrme models lead to the harmonic map ansatz J3], in 
which the original rational map of Houghton et al l* was 
replaced by the projector of the more general harmonic 
map of the CP^^^ model and then to the realisation 



that one could use more of such projectors, which repre- 
sent harmonic maps, to obtain further solutions (radial 
cases) or field configurations of the SU{N) models. Since 
then, these ideas have been used to other models (i.e. 
monopoles) or similar ones coupled to gravity. 

All these studies have relied on one fundamental as- 
sumption: the separation of the angular degrees of free- 
dom given in terms of a projector or a series of projectors 
and the radial dependence built in through the profile 
function or functions (for more projectors). 

However, it is easy to check that the profile function 
does not have to depend only on r; it could also depend 
on the angular variables. As long as the projectors repre- 
sent harmonic maps - the profile functions can be more 
general, and therefore provide a better approximation to 
the solutions of the model. This observation constitutes 
the essence of our improved harmonic map ansatz; 
the profile functions in addition to their dependence on 
r do depend, also, on z and z which allows for a more 
general angular dependence of the fields and density func- 
tions. 

In this paper we present our ansatz and reinvestigate 
the lowest baryon number field configurations of the 
SU{2) model. In particular we look at i? = 2, 3 and 4. 
In the subsequent papers we will look at larger values of 
B; the SU{N) models for N > 2; and the multiprojector 
case. We will also look at the gravitating cases involving 
skyrmions and monopoles. 

Tie Harmonic Map Ansatz and its Improved Version 

The SU{N) Skyrme action to be considered is given 
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by 



^tr (K, Kn + ^tr {[K„ K,] [K^\ K'^]) 



(1) 



where K,, = df^UU'^ for /x = 0, 1, 2, 3; U is the SU{N) 
chiral field and k, e are coupling constants. Although, 
we present the discussion for the general SU{N) model 
in this paper we will give the details of our studies for 
the SU{2) case. 

It is well known that, in order for the finite-energy 
configurations to exist, the Skyrme field has to go to a 
constant matrix at spatial infinity: U —* /as \x'^\ — > 
oo. This eff'ectively compactifies the three-dimensional 
Euclidean space into S^ and hence implies that the field 
configurations of the Skyrme model can be considered as 
maps from S^ into SU{N). 

This compactification leads to the existence of a con- 
served topological current yielding the topological charge 
to be identified with the baryon number B defined as 
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and gA'i'a/a jg ^]^g (constant) fully antisymmetric tensor. 

The starting point of the further discussion is the in- 
troduction of the coordinates r, z, z on R . In terms of 
the usual spherical coordinates r, 0, </> the Riemann sphere 
variable z is defined by: z = e"^ tan(6'/2). Then, the 
Skyrme action becomes [l| 
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while the baryon number takes the form 

B = / tr {Kr [K„ Kg])dr dz dz. 



(5) 



Since we are interested in the static field configurations 
in what follows we assume no t dependence. 

Next we consider the harmonic map ansatz which in- 
volves assuming that the Skyrme field is of the form [3| 



where P is a N x N hcrmitian projector which depends 
only on the angular variables (z, z) and h is the profile 
function which depends, at least, on r. Note that, the 
matrix P can be thought of as describing a mapping from 
52 into CP^-^ defined as 
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where V^ is a iV-component complex vector (depending 

on z and z). For N = 2 and V = (1, ,f{z)) where f{z) is 



a rational function we recover the rational map ansatz of 
Houghton et al [J. 

Following J3] , we define a new operator P+ by its action 



N 



on any vector u G (C as 



p+v = d. 



V {v'^ dzv) 



Note that, when V — V{z) (only) then 

P+V^V"! 



and so 
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where the subscripts denote partial derivatives. 

For © to be well-defined at the origin, the profile func- 
tion h, as a function of r, has to satisfy h{Q) = ir while the 
boundary value [/—*■/ at r = 00 requires that h{oo) — 0. 
As shown in [l|, Q{ , an attractive feature of lO is that it 
leads to a simple expression for the energy density which 
can be successively minimized with respect to the pa- 
rameters of the projector P and then with respect to 
the shape of the profile function h. This procedure for 
h = h{r) gives good approximations to multiskyrmion 
field configurations 0,|2|- 

In what follows, this ansatz is "improved" by allowing 
the profile function h to depend on z and z in addition to 
its r dependence. The new ansatz is consistent with the 
partial factorisation of the field in the sense that it still re- 
duces the problem to having to solve one equation for one 
function - namely h. Had we taken an ansatz in which, 
say, the parameters of the projector P depended on r, 
such a simplification would not have taken place. Thus, 
this modification is "nontrivial" and that is the reason 
for calling it an improved harmonic map ansatz. Note 
that, h has to be real implying that h — h{r, jzp, 4jp); 

while, at the origin, we require that h{0, jzp, n^) = tt. 

The action Q, due to ^ ioi h = h{r,z,z) and us- 
ing the aforementioned properties of the harmonic maps, 
becomes 
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where 
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FIG. 1: The one skyrmion profile function showing no 6 de- 
pendence (left) and the two skyrmion profile function ob- 
tained from the improved harmonic and rational map ansatz 
(right). 



while the baryon number @ coincides with the expres- 
sion for the topological charge of the CP'^~^ sigma model 
(up to an overall profile dependent factor) since 
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It is easy to see that /p°° sin^ hh^dr = h{r = 0)/2 = 7r/2. 
Next, for convenience of our numerical simulations, we 
use spherical coordinates (r, 9, (p) and introduce the di- 
mensionless coordinate x = enr. Variation of (|ll|l with 
respect to h gives the equation of motion 
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where G 
only. 
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function of 9 and 



SU{2) B = l,...,4Bar7ons 



First, to test our approach, we have calculated the pro- 
file function h for one skyrmion (i.e. for B = 1) and found 
no or dependence (as expected), as FigureHindicates. 
The total energy is 1.232. 

Axially symmetric skyrmions with baryon number B > 
1 can be obtained from vectors of the form [l| 
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(15) 



In this case, the action (|ll|l depends explicitly on 9 but 
not on (f>. Consequently, we find 9 dependence of the 
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FIG. 2: The skyrmion profile function h{r, 6) with axial sym- 
metry obtained from the improved harmonic and rational map 
ansatze for B = 3 (left) and B = 4 (right). 
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FIG. 3: The B = 3 skyrmion profile function h{r = 1,0, 0) 
with tetrahedral symmetry obtained from the improved har- 
monic map ansatze. 



corresponding profile functions h, as shown in Figures ^ 
and 121 The energies are given in Table 1 and compared 
with those obtained from the rational map ansatz. For 
B = 2 the energy is 2.382 which is closer to the value 
2.342 of the "exact" solution (obtained numerically by 
solving the full equations |J]) and lower than the value 
2.416 obtained by the rational map ansatz [l|. [In this 
case, h is 4> independent (as expected)]. For B > 2 the 
solutions (of this class) correspond to saddle points of the 
energy. Obviously, the improved harmonic map ansatz 
yields lower energies than the rational map ansatz, as 
can be seen in Table 1. 

Next, we consider skyrmions with platonic symmetry. 
For B — 3 and i? = 4 we let the vector V to be 



(/) V = (y3iz^-l,z{z^-iVi)) 
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All these expressions come from [ij - who established 
their form by minimizing the projector part of the ac- 
tion. This time, the corresponding action Hll|l depends 
explicitly on 9 and 0. 

First we consider the question whether we can obtain 
an improvement of the energy by considering h to depend 
(only) on 9, in addition to r. We have performed the in- 
tegration of the action over the azimuthal angle </> which 
has given us an effective action, which still contains an 
explicit 9 dependence. The variation of this effective ac- 
tion with respect to h{r, 9) led to a partial differential 
equation for it. We have found that the 9 dependence 
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FIG. 4: Same as Fig. I^for 5 = 4 with octahedral symmetry. 



of h is very small for B — 3 while for i? = 4 it is more 
pronounced, but still small compared to the axially sym- 
metric solution. As shown in Table 1, the improvement of 



J 



the energy is rather small, compared to the improvement 
of the energy for the axisymmetric solutions. 



This result suggests that a considerable improvement 
of the energy can only be achieved if we allow h to depend 
also on (/), in addition to r and 0. The variational equation 
then yields a second order partial differential equation 
involving the independent variables x,6, cf). We solved 
this equation for B = 3 and B = A and found indeed a 
further improvement of the energies (see Table 1). For 
B = 3, respectively, B = A the deviation from the exact 
value is only « 2 %, respectively, « 1 %. To demonstrade 
the angular dependence of the profile function we plot in 
Figures|niand^/i(a:, 9, (f>) for fixed x = enr = 1 for i? = 3 
and S = 4, respectively. 
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TABLE 1: Energies per skyrmion (i.e. E/B) of multiskyrmion configurations obtained by the rational and the 
improved harmonic map ansatz (with and without (f> dependence) for B — 1, . . . , 4 in comparison with the energies of 
the "exact" solutions. 
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Conclusions 

In this letter we have pointed out that the rational map 
ansatz of Houghton et al [l| can be further improved. The 
improvement involves the allowance of the profile func- 
tion h to depend on z and z in addition to r. Of course, 
the rational map ansatz is already a very good approx- 
imation; hence our improvement is only modest but, as 
we have shown (in the cases studied) the improvement is 
nonnegligible. To get our values we first have assumed 
only 9 dependence of h; however in order to get better 



values for the Skyrmions with platonic symmetries we 
went further and allowed, also, its cj) dependence. Our 
modification of the ansatz is not restricted to one projec- 
tor. It is easy to see that if we want to obtain low energy 
configurations of the SU{N) model we could use more 
projectors 2] (up to iV — 1 for the SU{N) case) with the 
corresponding profile functions dependent on r, 9 and 4>. 
We are also looking at such further applications of our 
ideas: to systems involving SU{N) Skyrme models for 
N > 2 and to systems involving gravitational fields. 
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